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ABSTRACT 

We present optical and infrared spectroscopy of V1309Sco, an object that erupted in 2008 in a stellar-merger event. During the 
outburst, V1309 Sco displayed characteristics typical of red transients, a class of objects similar to V838 Mon. Our observations were 
obtained in 2009 and 2012, i.e. months and years after the eruption of V1309 Sco, and illustrate severe changes in the remnant, mainly 
in its circumstellar surroundings. In addition to atomic gas observed in earlier epochs, we identified molecular bands of TiO, VO, H 2 O, 
ScO, AlO, and CrO. The infrared bands of CrO we analyse are the first astronomical identification of the features. Over the whole 
period covered by our data, the remnant was associated with a cool (<1000 K) outflow with a terminal velocity of about 200kms“*. 
Signatures of warmer atomic gas, likely to be still dissipating the energy of the 2008 outburst, dramatically decreased their brightness 
between 2009 and 2012. In addition, the source of optical continuum disappeared sometime before 2012, likely owing to the formation 
of new dust. The final stage of V1309 Sco’s evolution captured by our spectra is an object remarkably similar to an older red transient, 
V4332 Sgr. In addition to providing a detailed view on the settling of the eruptive object, the observations presented here reinforce 
the conclusion that all the Galactic red transients are a manifestation of the same phenomenon, i.e. a stellar merger. The late spectra 
of VI309 Sco also suggest peculiarities in the chemical composition of the remnant, which still need to be explored. 
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1. Introduction 


VI309 Sco was discovered in outburst in 2008 by K. Nishiyama 
and F. Kabashima in Japa n and independently by Guoyou Sun 
and Xing Gao in China (iNakano et al.l l20()^ : it was initially 
classified as a nova (Nova Scorpii 2008). Detailed spectroscopic 
observations of the object during the outburst and in its early 
decline showed, h owever, that it did not follow the behaviour 
of classical novae (iMason et alJl201(ih . Low outflow velocities 
of a few hundred km s ', a presence of low-excitation gas, and 
- most importantly - evolution towards cool phase character¬ 
ized by an M-type continuum, pointed towards classification of 
VI309 Sco as a ‘red nova’ or a red transient (the group is also 
known as intermediate-luminosity optical transients or merge- 
bursts). So far, the best studied object in this group is V838 Mon, 
which erupted in 2002, and therefore the red transients are also 
known as V838-Mon-type objects. This newly recognized class 
of variables has been considered to eru pt as a result of a stellar 
collis i on and consequent coalescence (ISoker & Tvlen^ I2003L 


l2006t iTvlenda & Sokeill20()^ ISoker & Tvlen^2007 ). The case 

of VI309 Sco provided the best and direct observational ev¬ 
idence that red transients indee d erupt as a result of a stel¬ 
lar merger. iTvlenda et al.l (l201 ih presented photometry of the 
object going back nearly seven years before its outburst and 
showed that it was an eclipsing binary system with a period 
of about 1.4 days. This period was shortening at an exponen¬ 
tial rate before the 2008 eruption showing directly, for the first 
time, a spiralling-in process that eventually led to a merger. 
Given that VI309 Sco showed observational characteristics typ¬ 


ical of red transients and the strong evidence for the merger in 
2008, the object is considered the ‘Rosetta stone’ linking the 
class of eruptive variables with mergers, solely on strong ob¬ 
servational grounds. Five galactic objects belong to th is group, 
including OGLE-2002-BLG-360 dTvlenda et alJl2(J13h: the old¬ 
est of them are V4332 Sgr, which erupted in 1994 (iMartini et all 
1 19991 : ITvlenda et al.l l2005l). and possibly CK Vul with a very 
ancie nt eruption in 1670 dTvlenda et al.l 120131 : iKamihski et al.l 
l2015h . Many of the recognized red transients have been observed 
in ne arby galaxies, e.g. M31RV dRich et ^ 19891: Mould et al 


1990h. M31N 2015-Ola dKurtenkov et alJl2015nWilliams et al 


2015), and M85 OT2006-1 dKulkarni et alj2007 ). Red transients 


(hence mergers) are now considered to be a more common phe¬ 
nomenon than thought a few years ago, with a n average Galactic 
rate of 1 event per 2 yr dKochanek et alJl2014l) . Many will be ob¬ 
served with the current and future sky surveys. 


Observing red novae is the best way to study the 
merge r process, but theo retical effort has also been under¬ 
taken dNandez et alJ 1201^ . After their outbursts all red tran¬ 
sients develop complex circumstellar environments abundant in 
molecules and dust. Their characteristics are linked to the crucial 
information on the physics of a merger, for instance the mass 
and angular momentum lost from the system during the violent 
event. Here, we present spectra obtained with Xshooter on the 
Very Large Telescope (VLT) demonstrating spectral evolution 
of V1309Sco in 2009 and 2012, i.e. months to years after its 
spectacular outburst. The spectra cover the range between about 
300 nm and 2.48/tm providing a very comprehensive view on 
the post-merger remnant of VI309 Sco. This study is a continu- 
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ationof spectroscopic monitoring of this object bv iMason et alJ 
(1201 Ol) . whose last published observations were obtained on 21 
April 2009, i.e. 227 days after the maximum light. 

2. Observations 

All observations of V1309Sco were conducted with Xshooter 
(IVernet et al.ll201 iT) . which is a spectrograph mounted on Unit 
Telescope 2 of VLT. Xshooter registers spectra in three frag¬ 
ments simultaneously: ultraviolet within 300-560 nm, visual 
within 550-1020nm, and near-infrared between 1020 and 
2480 nm. We refer to the three spectral ranges as UVB, VIS, 
and NIR, respectively. 

All the spectra were acquired in the nodding mode to bet¬ 
ter correct the NIR data for the sky background. The length of 
the nodding throw was 5". The position angle of the nodding 
direction was different for different dates. The slit length in all 
the observations was 11". The observations were conducted with 
the atmospheric dispersion corrector inserted in the optical path. 
The detector was always used in the high-gain mode. A technical 
summary of the observations is provided in Table[T] 

Commissioning (data The hrst observations were obtained dur¬ 
ing one of commissioning runs of Xshootefl The observations 
were obtained on 4 May 2009 and were accompanied by see¬ 
ing of 077. In addition to V1309Sco, the slit covered a nearby 
star ~2" north from the eruptive object. The number of expo¬ 
sures varied for each arm (see Table[T]) and the exposure times 
between 200 and 900 s resulted in the total time on source of 
about 50 min. Even the shortest exposures in the VIS arm of 
400 s resulted in the saturation of the Her line. We used the slit 
widths of 1" for the UVB arm and of 079 for the two other arms. 
This provided us with spectra at the nominal resolving power of 
R = /t/A/l=5100, 8800, and 5300, respectively. 

The telluric standard star HIP89486 (spectral type B9 V) was 
observed immediately after the science observations, while the 
spectrophotometric standard star EG274 was observed a few 
hours later at the end of the commissioning run. 

Science-verification data Later in 2009, additional observa¬ 
tions were obtained within project 60.A-9445(A) (PI: E. Ma¬ 
son), which was part of the science veriheation of Xshooter. 
VI309 Sco was observed twice on 12 August and 29 September. 
Seeing during the observations was of 170. The slit widths and 
resolution reached were the same as in the earlier commission¬ 
ing data, but the detectors of the UVB and VIS arms were binned 
in the spatial direction (1x2 pix). The slit was aligned with the 
parallactic angle at the beginning of the exposure. A series of ex¬ 
posures was taken (Table[T]i adding to the total integration time 
of about 33 min both in August and September. 

The telluric standard HIP088857 (spectral type B3V) was 
observed at the two science-veriheation runs. The spectrophoto¬ 
metric standard star BD-i-17 4708 (spectral type sdE8), was ob¬ 
served only during the August run, but we used it to calibrate 
data from both runs. 

The 2012 data Another set of observations was obtained on 
22 April 2012 within a service-mode program 089.D-0041(A) 
(PI: T. Kaminski). Two spectral setups were used. The hrst em¬ 
ployed the widest available slits of 5" to minimize slit losses 

* ESO programme 60.A-9024(A), see 

https://WWW.eso.org/sci/activities/vltcoimn/xshooter.html 


and to provide a spectrum suitable for accurate hux calibration. 
The spectral resolution was between 1900 and 3200 in the three 
spectral ranges. The second setup was aimed to provide spectra 
at a higher resolution and a set of narrower slits was used: 170, 
077, 074 for UVB, VIS, and NIR, respectively. This provided us 
with spectra at nominal resolutions of R-4 300 for UVB, and of 
11000 for VIS and NIR. The average seeing during the 2012 ob¬ 
servations was of 075. The 2012 observations were considerably 
shorter than those in 2009 (cf. Table[T]i. 

During the 2012 observing run, the telluric standard 
HIP067796 (spectral type B2II) and the spectrophotometric 
standard EG274 (spectral type DA) were observed on the same 
night as VI309 Sco. A telluric standard was observed with a slit 
width of 079 in VIS and NIR arms, and the Hux standard was 
observed with the wide 5"slits. 


3. Data reduction 

The commissioning and science-verihcation data from 2009, as 
test observations, required special care when basic image pro¬ 
cessing was applied. The primary data reduction was obtained 
using the esorex ESO pipeline (version 1.2.2) and included Hat- 
Held correction, echelle-order rectiHcation, and wavelength cal¬ 
ibrations. Eor the commissioning data, all the necessary calibra¬ 
tion frames were taken the day following the science observa¬ 
tions. However, spectra of calibration lamps were underexposed 
and could not be used alone. We therefore used corresponding 
frames taken a few months later. We found the same strong emis¬ 
sion lines in the underexposed and well-exposed spectra of cal¬ 
ibration lamps; based on their comparison we can exclude sys¬ 
tematic errors in the wavelength calibrati on due to possible (bu t 
unlikely) changes in the spectral format (iMartavan et alj|^14l) . 

Eurther processing of all the 2009 data — including the order 
extraction, order merging, telluric correction, and flux calibra¬ 
tion — was performed using standard procedures in IRAE. We 
applied the telluric correction in IRAE only to the NIR-arm spec¬ 
tra. The telluric correction in the VIS re gion was applied t o 
the flux-calibrated spectra using molecfit (ISmette et alJl2015T) . 

This tool corrects for the telluric absorption by performing a 
full radiative-transfer simulation of the atmospheric absorption 
bands at weather conditions specifle for the science observations. 

The telluric correction was calibrated on the visual O 2 bands, 
which do not have any circumstellar components. The contribu¬ 
tion of telluric water vapour was then calculated using a standard 
H20-to-02 abundance ratio. The applied procedure left residua 
below ~10% of the original telluric absorption. 

The basic data reduction of the 2012 data (including flat- 
field corrections, wavelength, and flux calibrations) were ob¬ 
tained usi^ the standard Xshooter pipeline under the Reflex en¬ 
vironment However, owing to a severe contamination by Held 
stars, the spectral extraction had to be performed manually and 
separately for each of the nodding phases. This was done in 
IRAE. Eor reference, the 2012 data were fully reduced using 
the pipeline in different versions, and we also made use of the 
Xshooter Data Release especially for the presentation of the 
NIR-arm data. The cumbersome spectral extraction owing to the 
Held crowding and an unfortunate position angle of the nodding 
direction practically exclude any reliable flux calibration of the 
2012 data. Therefore, the data obtained with the wide slit are 
not presented here. We corrected the 2012 spectra for telluric 
absorption using molecfit. Because the NIR spectrum has no 

^ http://WWW.eso.org/sci/software/pipelines/ 

^ http: //WWW. eso.org/ohserving/dfo/guality/PHOENIX/XSHOOTER/pr 
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Table 1. Technical details of the observations. 


UT date 

Days 

since 

max. 

Slit 

PA'* 

(deg) 

exposure 
time (s) 

UVB arm 
slit 

width (") 

binning 

(pix) 

VIS arm 
exposure slit 
time (s) (") 

binning 

(pix) 

NIR arm 
exposure slit 

time (s) (") 

binning 

(pix) 

2009-May-04 

240 

1.7 

600x2 

1.0 

1x1 

600x2 

0.9 

1x1 

200x6 

0.9 

1x1 

2009-May-04 

240 

1.7 

900x2 

1.0 

1x1 

400x4 

0.9 

1x1 

450x4 

0.9 

1x1 

2009-Aug-13 

341 

51 

1000x2 

1.0 

1x2 

1000x2 

0.9 

1x2 

(300x3*)x2 

0.9 

1x1 

2009-Sep-29 

388 

51 

1000x2 

1.0 

1x2 

1000x2 

0.9 

1x2 

(300x3*)x2 

0.9 

1x1 

2012-Apr-22 

1324 

-72 

200 

5.0 

1x1 

200 

5.0 

1x1 

268 

5.0 

1x1 

2012-Apr-22 

1324 

-72 

394x2 

1.0 

1x1 

394x2 

0.7 

1x1 

461x2 

0.4 

1x1 


Notes. Position an gle measured from North through East. *** We did not obtain six separate exposures but rather set NDIT parameter 
dMartavan et alJl?0l4) to 3, so that single frame was stored for each nodding position and corresponds to a pre-processor average of three in¬ 
tegrations of 300 s each. 


significant continuum, the molecf it correction procedure could 
not be applied for it. 

The wavelength calibration of the science-verification and 
the 2012 data is accurate within an rms of 0.002 nm in the 
UVB and VIS range, and 0.01-0.02 nm for the NIR arm. All 
spectra and radial velocities are presented here in the heliocen¬ 
tric rest frame, Vf,. The correction from Vk to the local stan¬ 
dard of rest (LSR) velocity at the position of VI309 Sco is 
Vlsr=V/,- 1-9.775 km s^'. 

In Fig. [T]we show the Xshooter acquisition image illustrat¬ 
ing the location of the slits in the 2012 observing run. The image 
was obtained in the V filter in 2012. The wide slit encompassed 
several bright stars in the field, in particular the star 2" north 
from the position of VI309 Sco. This source dominates or is 
comparable to continuum brightness of VI309 Sco from optical 
to NIR wavelengths in the time period covered by our observa¬ 
tions (Kaminski e t al. in prep.). Images o btained with the Hubble 
Space Telescope (iMcCollum et al.ll2014l) also resolve the source 
seen in the Xshooter image at the position of VI309 Sco into 
at least three separate sources, so also the observations with the 
narrow slits gave spectra considerably affected by the field stars. 
An inspection of our earliest-epoch data, i.e. when VI309 Sco 
was brightest, show that even then the spectra were contaminated 
by the field stars. Because of this blending effect and the rough 
absolute-flux calibration, we focus here on an analysis of spec¬ 
tral features in the spectrum of VI309 Sco, not the continuum. In 
particular, we subtracted the false continuum of the 2012 spec¬ 
trum as a low-order polynomial. We did not correct the spectra 
for interstellar extinction. 


4. Identification and analysis of spectral features 

The VIS and NIR spectra with the proposed identification are 
shown in Figs.[3-(6] (The UVB spectrum will be presented 
elsewhere.) To identify atomic emission lines, we used the 
data base of NIST (Nation al Institute of Standards and Technol¬ 
ogy ;|&a^daet^02W3)- The location and s tructure of mol ecu- 
lar bands was simulated mainly in pgopher (IWesternll2014ll us¬ 
ing spectroscopic constants and literature data described in detail 
in AppendixlAl 

As shown in Fig.|3 the spectra were obtained during the 
fading part of V1309Sco’s light curve and the continuum flux 
changed by a factor of about 24 from the earliest to the latest 
spectrum. We analyse the spectra emphasizing the time variabil¬ 
ity of spectral features in that period. 



Fig. 1. Location of the VIS slits in the 2012 observations. 



Fig. 2. Dates of Xshooter observations marked on the light curve o f 
VI309 Sco as observed by OGLE in the / filter dTvlenda et alll20lih . 
The light curve before and after MJD=5000 is shown in different scales 
as indicated on the left (blue) and right (red) axes. 

4.1. UVB range 

Only very weak continuum rising towards longer wavelengths 
was detected in the UVB range. Furthermore, weak absorption 
features are present longward of about 490 nm, with the Mg i 
triplet (Mgb, 516.7-518.4nm) and H/I being the most promi¬ 
nent examples (the former seen clearly only in the 2012 spec- 
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tmm). The continuum and absorption spectrum most likely be¬ 
long to the field star(s) of spectral types close to G-K. All the 
UVB spectra are very rich in emission lines, which undoubtedly 
do belong to V1309 Sco. 

Each spectrum acquired in 2009 has over 60 spectral fea¬ 
tures. Although obtained weeks and months apart during the 
very steep dimming part of VI309 Sco’s light curve (Fig.|2]), the 
2009 spectra appear very similar to each other when scaled to 
the same intensity. No spectral features have been found above 
the noise levels shortward from about 380 nm, except for the ear¬ 
liest Xshooter spectrum where emission lines are seen down to 
at least 320 nm. 

In the 2012 spectrum, which has the lowest signal-to-noise 
ratio, only the strongest of the lines seen in 2009 are clearly vis¬ 
ible and many of them are narrower in 2012 (all UVB spectra 
have a very similar resolution). A few new features are present 
in the latest spectrum, e.g. near 480.5 and 490.6 nm. 

In all the spectra, the most prominent emission lines are the 
Balmer features (except for the 2012 spectrum where they are 
absent), individual lines of [O i] and Mg i], and two lines of S ii. 
In addition, several multiplets of Fe i and Fe ii (allowed and for¬ 
bidden) and the H&K lines of Ca ii were securely identified in 
the spectra. However, more than half of the observed spectral 
features could not be unambiguously identified. Many of those 
unidentified features are at wavelengths of known lines of Fe, Ti, 
Cr, V, Ni, and their ions, but not all lines in the given multiplet 
could be traced in the spectrum, hampering a definite assign¬ 
ment. Judging from the relative intensities of the identified lines 
of Fe I and Fe ii, the emission spectrum does not arise from is- 
tothermal gas. Many lines may be produced by scattering which 
results in line ratios strongly deviating from those expected from 
emission of gas in thermal equilibrium. 

When the emission spectra from the different epochs are 
compared in absolute scale, a clear decrease in line fluxes 
with time is apparent. The continuum-free integrated flux in 
the 385-559nm range dropped from 13.9 (May 2009) to 2.7 
(August 2009) and then to 1.9 (September 2009) in units 
of 10 ergcm^^ s '. The same measurement for the 2012 
spectrum, although with very uncertain flux calibration, gives 
6.3xl0^*^ergcm“^s '. The decrease in the average line fluxes 
is therefore steeper than the corresponding drop in /-band con¬ 
tinuum documented by the OGFE light curve (Fig.|2]i. The dif¬ 
ferent lines fade in a slightly different rate. In particular, the Hi 
lines fade much faster than those of all the observed metals; for 
Hy6, the measured flux in units of 10 ergcm^^ s ' goes from 
39.3 (May) to 2.7 (August) and to 2.2 (September) in 2009 spec¬ 
tra, and the line is absent in the 2012 spectrum. The Can H&K 
lines faded very quickly as well. 

Although spectral features are present close to the predicted 
positions of heads of the B - X band of AlO (e.g. near 485 nm 
and 507.5 nm) and the a (0,0) band of TiO (at 516.5nm), there 
is no convincing evidence for a presence of molecular emission 
bands in the UVB range. 

4.2. VIS range 
4.2.1. Atomic lines 

In the VIS range of the four Xshooter spectra, we found only 
a few atomic lines, among which the strongest are the Ha line; 
lines of [Oi] (630.0 & 636.4nm) and [Fen] (e.g. 544.0, 564.0, 
569.6 nm), the resonance doublets of Ki and Rbi; the semi- 
forbidden line of Cai (seen outside the Ha line only in 2012); 


and the resonance pair of Cs i lines (684.9 & 689.5 nm), absent 
in 2012. 

While the 2009 spectra are undoubtedly dominated by the 
emission and absorption features of VI309 Sco itself, the spec¬ 
trum obtained in 2012 - when the object was substantially 
weaker in continuum (~20 times in I) - is contaminated by spec¬ 
tra of the blending objects that have a comparable brightness to 
V1309SCO (Fig.[B. The 2012 spectrum contains narrow absorp¬ 
tion features of which the Nai doublet at 589 nm and the Can 
triplet near 860 nm are most clearly recognizable. These absorp¬ 
tion features most likely come from the blending stars whose 
spectral type must be close to G-K, considering the strengths 
of the mentioned lines. At a closer inspection, the Can triplet 
also appears to have weak and blue-shifted (with respect to the 
absorption line centre) emission component, which can be as¬ 
cribed to VI309 Sco. This is supported by the presence of pure 
emission lines of Ca ii in the earlier epochs. 

The time evolution of most emission lines is linked to what 
was described above for lines in the UVB range - most of them 
quickly fade with time. These include Ha and features of Fe ii, 
Ca II, and [O i]. This group of lines gets fainter with time in terms 
of absolute flux units and also with respect to the continuum. 

An interesting behaviour is observed in resonance lines of 
alkali metals. The line of Csi near 894.3 nm gets weaker with 
time over 2009 to completely disappear in 2012. The /17800 line 
of Rbi in 2009 was getting fainter with time in absolute flux 
units, but was gradually gaining its strength with respect to the 
local continuum (Fig.O. In the 2012 spectrum, which however 
has very poor flux calibration, the 47800 line is comparable or 
even higher than in the last spectrum from 2009. The other line 
of the Rbi doublet, 47944, although clearly present in 2012, is 
not seen at all in the 2009 spectra, possibly because it is ab¬ 
sorbed by the saturated absorption band of VO. This implies, 
that the source of atomic emission must be internal to (behind) 
the cloud of absorbing molecular gas. The optical doublet of Ki 
shows very similar behaviour to the 47800 line of Rb i - over the 
2009 observations, it fades on a pace lower than the decay of the 
underlying continuum but in 2012 the doublet is clearly stronger 
than in any of the 2009 spectra. 

The optical doublet of Nai shows a different behaviour. In 
2009 spectra, the lines show weak net emission despite being 
absor bed by interstellar absorption components (cf. lMason et al.l 
1201 oh and photospheric lines of the blending field stars. The 
emission gradually fades in time during 2009 observations. Only 
very weak emission is seen in the 2012 spectrum, comparable to 
that observed in the first 2009 spectrum. Because the Na i usually 
produces stronger lines than Ki in a wide range of astrophysical 
environments and physical conditions, the observed characteris¬ 
tics of those two doublets in V1309 Sco is rather peculiar. A sim¬ 
ilar peculiarity has been noticed in V4332 Sgr and is explained 
as the r esult of the scattering of light of a source that is hidden 
by dust (iKaminski et alJl^lOUTvlenda et al.ll2014h . This feature 
bears crucial information on the emission mechanism and source 
geometry as discussed in more detail in Sect.|6] 

4.2.2. Molecular bands and pseudo-continuum 

The VIS range is abundant in molecular bands. The 2009 spec¬ 
tra are very similar to each other and show a cool continuum 
absorbed in bands of simple oxides which are typical for late M- 
type stellar photospheres. These bands belong mainly to the elec¬ 
tronic systems of TiO (y, y', e, 6) and VO (B - X) (Fig.[3]l. Ad¬ 
ditionally, above 890 nm, we observe strong absorption of water , 
most clearly in the pm ro-vibrational band (iTanaka et al.lll982h 
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Fig. 3. Identification of spectral features in the VIS range. The spectra from May 2009, August 2009, September 2009, and April 2012 are shown 
with magenta, green, blue, and red lines, respectively. The spectra were corrected for telluric absorption, but location and structure of telluric bands 
is shown at the bottom part of each panel (grey line). A continuum was subtracted from the 2012 spectrum. All spectra were smoothed and are not 
in scale with respect to each other. The top parts of the profiles of the Her and the [O i] lines in the 2009 spectra were clipped for clarity and are 
marked with the scissors (some profiles are compared separately in Fig.|7j. Simulations of the (0,0) band of the VO B- X and CrO A' - X systems 
are shown in the lower panel below the 2012 spectrum (in cyan and black, respectively). They were generated with pgopher at an excitation 
temperature of 300 K (see AnnendixlAt. 


whose upper level is 7650 K above the ground. The weaker vi¬ 
sual bands of water, a (~723nm) and // (~822nm), are also 
present but less conspicuous. We are confident the absorption in 
the pcTT bands is intrinsic to V1309 Sco and not to the Earth’s at¬ 
mosphere because the lines are significantly broader than the in¬ 
strumental profile and are shifted to Vh of —160+20km s“'. Ad¬ 
ditionally, we compared the observed spectrum (corrected for the 
telluric features) to a grid of simulations of the pen band at dif¬ 
ferent excitation temperatures. The observed shape of the head 
near 927 nm can only be formed in excitation temperatures sig¬ 
nificantly higher than those in the Earth’s atmosphere, as illus¬ 
trated in Eig.|4] At temperatures <380 K strong abortion is seen 
only above 932 nm, i.e. 6 nm longward from the observed head. 
The location of the band-head can be very well reproduced for 
850 K> Tex >1000K; this can be treated as the first rough con¬ 
straint on the temperature of the absorbing circumstellar vapour. 
There are other opacity sources overlapping with the water band, 
but - owing to the dominant contribution of water - they could 
not be conclusively identified. In particular, emission-like fea¬ 
tures at 966.9, 976.4, 989.9, and 1004.6 nm may be related to 
the vibrational band of OH (Eig.|5]l. 


With respect to the continuum, the absorption bands were 
deepest in the earliest observation of May 2009 and were con¬ 
siderably weakened in the August/September 2009 observations. 
In most of the bands, we observe consequent weakening of the 
bands depth with time. Despite the imperfect flux calibration, we 
are confident that this effect is real and not related to an imperfect 
correction for the instrument response function. Eor instance, the 
shape of the absorption through of the VO B — X (0,0j3 band 
(in 760-810nm) changes with time as expected for consequent 
lowering of the column density of absorbing gas — there are al¬ 
most no changes near the band-head while significant changes 
are seen in the tail of the band where opacity is lowest. 

We made an attempt to derive the spectral type for the 2009 
observations assuming the pseudo-continuum is produced by a 
stationary stellar-like photosphere. By comparing the observa¬ 
tions to simulated spectra of th e MARCS stellar-atmosphere 
models (iGustafsson et al.l l2008h . however, we found that the 
shape of the pseudo-continuum strongly deviates from any M- 
type photosphere with a standard chemical composition (i.e. 


To assign the rovibrational bands, we use the convention (Uup,tJio„) 
where the numbers are the upper and lower vibrational states. 
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Fig. 4. Water absorption band in V1309 Sco (black) is compared to sim¬ 
ulations of the pcTT band at different excitation temperatures (see Ap- 
pendixlAt. The observed band cannot be formed as a telluric feature 
because the structure of the hand at air temperatures does not form a 
clear hand-head near 926 nm. All simulations were smoothed to fit the 
width of the observed features. Although water dominates, there are ad¬ 
ditional opacity sources in this wavelenght range, e.g. the TiO e system. 


close to solar). Most apparent in V1309Sco is the greatly en¬ 
hanced absorption within the B - X system of VO and the pen 
band of H 2 O. A similar effe ct in VO bands has bee n reported for 
the spectrum of V4332 Sgr (iKamihski et ahllToiOh . while strong 
optical abs orption of water was se en in another red transient 
V838 Mon (iLvnch et al.ll2004l2007h . Such strong absorption in¬ 
dicates the presence of cool circumstellar gas surrounding the 
stellar remnant. Although we compared the spectrum to a broad 
range of photospheric models applying different values of red¬ 
dening, we were not able to assign a spectral type in the classi¬ 
cal sense. Even if the spectrum arises in a stellar-like object, its 
composition may be very different from a normal star and spec - 
tral typing would then inevitably fail (cf. IKamihski et al]l 201 (]ll . 
Nevertheless, the overall slope of the observed spectrum is close 
to standard spectra generated for effective temperatures <3300 K 
or spectral types later than M5. Additionally, it is clear that the 
spectrum is formed in an extended atmosphere, i.e. of a giant 
or a supergiant, because absorption features of simple hydrates 
(e.g. FeH, MgH, CaH, CaOH, CrH) - usually present in spectra 
of dwarfs - are absent. 

Some molecular bands also showed weak emission in the 
2009 spectra. The B - X system of VO is dominated by absorp¬ 
tion, but within the absorption dips weak emission features of 
VO are also seen, most clearly in the (1,0) band near 740 nm 

(Fig.El). 

The molecular absorption bands disappeared almost com¬ 
pletely in 2012 , which is directly related to the dramatic (appar¬ 
ent) fading of the continuum source (Fig. la. Most of the molec¬ 
ular bands seen in absorption in 2009 turned into a pure emission 
spectrum in 2012 (excluding bands of water). In addition to the 
well-known electronic systems of TiO and VO, the 2012 emis¬ 
sion spectrum contains weak but clear bands of CrO belonging to 
the B-X system, most clearly (and unambiguously) in the (0,1) 
band near 640 nm. This band has been observed i n emission only 
in on e other astrophysical source, V4332Sgr (IKamihski et alJ 
I 2 OIOI) . We also found evidence for weak emission of CrO in an¬ 
other electronic system. A' — X, which overlaps with emission 
of the TiO e (0,0) band near 843 nm (Fig. [3]). Because the bands 
overlap and have complex structure, this identification is only 


tentative at the moment; it is supported, however, by the pres¬ 
ence of CrO emission in other electronic systems within the VIS 
and NIR spectra (see below). 

Other less common species seen in emission in 2012 are the 
A ^ 113/2 - X^IA bands of ScO and YO (within 590-615nm, see 
Fig-El), which are however weak. No features of circumstellar 
water vapour are seen in the 2012 spectrum above the residuals 
of the telluric correction. 

4.3. NIR range 

Atomic features are relatively rare in the NIR range, with only 
recombination lines of hydrogen (e.g. Paschen Rfi - 6, possibly 
Brackett Hy) having a reliable identification. A narrow emission 
feature at 1065.6 nm seen in the two spectra from 2009 may be a 
transition of Ti i or [Sc ii]. The NIR range is dominated by molec¬ 
ular features that are described in detailed below, separately for 
each of the identified species. 

4.3.1. AlO 

In the NIR spectra from 2009, emission bands of the A ^Yl-X 
system of AlO were easiest to identify owing to their characteris¬ 
tic double structure and triangular shape of each of the spin com¬ 
ponents (ni /2 and 113 / 2 ). These infrared bands were f irst recog¬ 
nized in spectra of V4332 Sgr in iBaneriee et al.l(l2003l) . To better 
characterize the location of all the AlO bands in VI309 Sco, we 
performed a simple simulation of the A - X system which is in¬ 
cluded in Figs.|5] and ^Particularly prominent in the covered 
range are the (4,0) banqj, seen in all spectra and located long- 
ward of about 1222 nm, and the ( 1 , 0 ) band, located longward 
of about 1645 nm; the emission spike near 2010nm, although 
affected by a (telluric) feature of CO 2 , is also dominated by the 
(1,1) band of AlO. Other weaker features are present in the spec¬ 
tra, as can be seen by comparing the observations to the simula¬ 
tion in Figs. |5] and | 6 ] 

From the band shapes and the relative intensities of the ni /2 
and 113/2 components, it is evident that optical thickness close to 
the band-heads is very high, i.e. t » 1 , especially in the 03/2 
components. Our simulation does not account for such strong 
saturation effects. The presence of bands involving high vibra¬ 
tional states (at least up to v-6) and relative intensities of the 
bands suggest a high vibrational temperature of 1000-3000 K, 
while the band shapes away from band-heads strongly suggest 
a rotational temperature of a few hundred K. These tempera¬ 
ture constraints are very uncertain, however, and should be veri¬ 
fied by a more realistic excitation model with a proper treatment 
of optical thickness effects and going beyond the assumption of 
isothermal gas. 

In our latest NIR spectrum from 2012, only the (1,0) and 
(2,0) bands of 03 / 2 , and (3,0) of ni /2 are undoubtedly present. 
The spectral region that was dominated in 2009 by the (4,0) band 
of AlO evolved until 2012 into a broad feature, which- as we 
discuss below- is dominated by CrO emission. 

4.3.2. VO 

Also present in the NIR spectra, most clearly in the data from 
2009, is the A ^Tt-Y system of VO. The (0,1) band gives rise 
to an emission feature between about 1165 nm and 1200nm that 

^ Because our simulation does not include pert urbations of the 
X v'=4 state by A ^If u=9 dLaunila & JonssonI 19941) . bands involving 
v'=4 are slightly shifted and incorrect in intensity. 
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is recognizable by its characteristic eight separate peaks (two 
features in each quartet component), as illustrated by a simple 
simulation of the system in Fig.|5] The band overlaps with a 
much weaker emission of AlO in band (6,1). The red shoulder of 
the (0,1) feature of VO may be strongly contaminated by the next 
band of the Ai;=l progression, i.e. the (1,2) band, which is also 
expected to have eight peaks. Owing to the lack of appropriate 
spectroscopic constants, we were not able to include this higher 
band into our simulation. Only very weak emission of the (0,0) 
band and perhaps also of (0,2) are present in the 2009 spectra. 

The VO bands are also responsible for absorption features, 
most clearly in the (0,0) band in the spectrum from May 2009 
(commissioning). Most other bands of VO can be interpreted as 
a combination of emission and absorption features. This also ex¬ 
plains why the (0,0) band emission is weaker than in (0,1) in 
all the 2009 spectra. It is apparent that each of the eight sub¬ 
heads of the (0,0) band seen in emission has a corresponding ab¬ 
sorption component blue-shifted by about 200 km s '. This self¬ 
absorption profile implies that photons produced by the warm 
molecular gas responsible for the emission passes through a 
colder layer of gas. The emission band (0,0) of VO is weak - 
yet still recognizable - in the 2012 spectrum. The time sequence 
of the NIR spectra of VO illustrates the same changes from pure 
absorption profiles to emission, similar to what is seen in the 
visual electronic system of VO. 

4.3.3. CrO 

Using a simulation of the CrO (0,0) band of the - V^Il 
system, implemented in pgopher by C. M. Western, we found 
a very satisfactory match of the 1220-1270nm feature seen in 
2012 to the simulated band shape at rotational temperatures close 
to about 300 K. The quintet nature of this electronic band is re¬ 
sponsible for the broad and very characteristic shape; we note 
that the observed profile may be still somewhat contaminated by 
the bands (4,0) and (5,1) of AlO. This band is remarkably strong 
and dominates the emission flux in the whole spectrum covered 
by Xshooter, which is surprising given that this is the first as- 
trophysical identification of this band of CrO. We extended the 
simulation of the CrO A-X system to include the corresponding 
(1,0) and the (0,1) bands. Although covered by our 2012 spec¬ 
trum, the bands fall within absorption bands of water and could 
not be unambiguously identified. The bands may be weak owing 
to a low gas vibrational temperature in 2012. 

To our knowledge, this is the first astrophysical identification 
of this infrared CrO band in an astronomical spectrum. 

4.3.4. CO 

The 2009 spectra above 2.3 pm show the characteristic rotational 
pattern of the rovibrational (2,0) band of CO (in the ground elec¬ 
tronic state X 'S). We interpret all the features above 2.3 pm as 
emission of the Av-2 sequence, although some contribution from 
underlying CO absorption cannot be excluded. The shape of the 
profile suggests the head of the (2,0) R branch is not formed 
(or is highly saturated), which implies a rotational temperature 
below about lOOOK. However, rotational components from vi¬ 
brational levels above u-2, i.e. from (3,1), (4,2), etc., are seen 
and suggest a much higher vibrational temperature of at least 
2000 K. At the same time, however, there is no evidence for the 
presence of higher CO overtones in the 2009 spectra. By cross- 
correlating the observed spectra against our simulation of the 
Av=2 sequence, we found the velocity of the two spectra did not 


change within 2.4 km s“' in the 47 days between their exposures 
and that the mean radial velocity was -117 + 7kms“' (in the 
heliocentric rest frame). There is no significant deceleration of 
the flow. 

In the 2012 spectrum, only weak absorption features can be 
seen at the expected positions of the main band-head of the Aii=2 
sequence of CO and the presence of CO features is tentative 
at best. If present, the absorption features may originate from 
a spectrum of one of the blending stars because the location of 
the absorption bands corresponds to the band-heads formed at 
high rotational temperatures, typically met in photospheres of 
spectral types later than ~G0. The lack of circumstellar CO fea¬ 
tures in the late spectrum of VI309 Sco is surprising given the 
presence of strong and clear emission of rare and exotic species. 

4.3.5. H2 

In the 2012 spectrum where we expected to see CO bands, we 
observe one emission feature which we identify as the Q{\) (i.e. 
7=l-l)t;=l-0 rovibrational transition of H2 (labora¬ 
tory wavelength 2406.6 nm). There are no other features that can 
be unambiguously identified as H2 transitions. If this identifi¬ 
cation is correct, all other i) = 1 - 0 lines belonging to the S 
branch and those in the Q branch from rotational levels 7 > 1 
must be at least ~6 times weaker than 2(1). In particular, we do 
not observe the 1 - 0 5(1) (i.e. 7 = 3-1) line at a laboratory 
wavelength of 2121.8 nm which is expected to be of comparable 
intensity to 2( 1) in excitation conditions typical for most known 
H2 emission sources, i.e. in thermal equilibri um conditions, o r 
fluorescent- or shock-induced excitation (e.g. ISternber3ll989ll . 
Perhaps, the line is absent owing to an excitation mechanism that 
is not efficient above 1 7=1 making other lines very weak. The 
H2 emission extends over a range from -350 to -10kms * and 
has a nearly rectangular profile. 

4.3.6. H2O 

Similar to what is seen in the optical region, the NIR spectra 
show circumstellar features of the rovibrational bands of water in 
strong absorption. In fact, the water bands dominate the overall 
appearance of the NIR spectra in 2009. The lines are broad and 
blueshifted leaving no doubt they are of circumstellar origin. A 
simple simulation of the (vi, v'2, V3)=(l,l,l) band near 1130nm 
- which (i) has the highest-energy above the ground (6120K) 
among all the observed NIR bands, (ii) is least contaminated by 
residuals of the telluric correction, and (iii) is least saturated - 
gave us an excitation (rotational) temperature of about lOOOK, 
i.e. consistent with that derived for the optical pen band. The 
absorption bands to lower-lying vibrational levels, i.e. (1,0,1) 
and (0,1,1), are saturated and not suitable for a simple excita¬ 
tion analysis as ours. As in the VIS spectrum, there is no sign of 
water features in the 2012 data, most likely owing to the lack of 
the continuum source. 

4.3.7. Unidentified features 

Numerous features observed in the NIR spectra remain uniden¬ 
tified. The vast majority of them are broad and display rotational 
patterns (i.e. a high-frequency “wave”) that strongly suggests 
that their origin is molecular. By performing simple simulations 
of rovibrational or electronic spectra of numerous oxides (MgO 
A'n - V'n, CaO A' 'n - V and TiO (f) and 6 systems), hy¬ 
drates (e.g. FeH -F'^A and HF V - X 'S+), and other 
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compounds (e.g. ZrS a^A-b' and CN we ex¬ 

cluded some systems of these species as carriers of the observed 
features. Unfortunately, some of the bumps seen in the spectra 
are due to imperfect data reduction and instrumental artefacts 
whose deforming influence is often very difficult to asses in the 
NIR part of our Xshooter spectra. 

5. Line profiles and velocities 

The kinematical changes that the remnant underwent during the 
period covered by our observations can be traced most clearly 
in the profiles of atomic features. There are two major groups of 
atomic lines that show distinct type of profiles. Example profiles 
are shown in Fig.|2] 

Forbidden emission lines of metals, e.g. Mgi], [Feii], [Oi], 
and the Balmer lines, show pure emission profiles over the 2009- 
2012 period. As described earlier, the lines weaken and become 
narrower with time. For instance, the 715.5 nm line of [Fe ii] had 
a FWHM of 130kms“' in May 2009, llOkms^' in Aug./Sept. 
2009, and 57 km s“^ in 2012; similarly, the width of the strongest 
line of [Oi] dropped from 131, through 118, to 87kms"*. The 
[Oi] and Hi lines also became more asymmetric with time. 
These lines are centred at around-180km s ' with a scatter and 
systematic differences between the different species of a few 
kms *. In this group, the Balmer lines are special. The 2009 
profiles of Ha contain broad Forentzian-like wings that can be 
traced as far as ~500kms * from the line centres, most likely 
produced as scattering on electrons. The Ha profile in 2012 does 
not show such wings but instead has an extra emission compo¬ 
nent at —12kms“*. 

Another group of characteristic lines consists of strong res¬ 
onance doublets of alkali metals, i.e. of Ki, Rb i, Cs i, which all 
showed P-Cyg type profiles before 2012. The absorption com¬ 
ponent of those profiles was most pronounced in the earliest 
observations and was becoming weaker with time; in 2012 the 
absorption component was weak or absent in some of the pro¬ 
files. The lines of Cs i diminished until they completely disap¬ 
peared in 2012. The broadest in this group are lines of Ki and 
their absorption extends down to at least y/,=-366kms *. The 
emission profiles peak near -95 kms * and at least some of the 
profiles have multiple components, for instance one peaking at 
about -12 km s'*. The turn-over point from emission to absorp¬ 
tion, which usually indicates the systemic velocity (the kinemat¬ 
ical centre of an outflow), is between -190 and -155 km s ', de¬ 
pending on the line and the somewhat arbitrary definition of lo¬ 
cal continuum level. Then the furthest-reaching absorption seen 
in KI indicates an outflow terminal velocity of about 200 km s '. 

Among the lines of alkali metals the profiles of the Nai 
optical doublet clearly stand out. The emission is very weak 
at all epochs. The 2012 spectrum has clear absorption features 
blueshifted with respect to the faint emission components. As 
noted earlier, these features may be a composition of the pho- 
tospheric and interstellar absorption features in the spectra of 
the blending stars. In earlier epochs, the interstellar component 
could have reduced the strength of the doublet somewhat, but the 
lines are intrinsically weak. 

The H2 line, which shows pure emission profile, extends over 
the whole velocity range covered by all the other lines and is 
clearly different from all the lines described so far. It may contain 
several discrete emission components blending into this single 
broad profile. 

Considering the overall kinematical appearance of all ob¬ 
served atomic features, the most characteristic velocity - which 
may well be the systemic velocity - is at about Vsys.h = 



-600 -500 -400 -300 -200 -100 0 100 200 300 

Heliocentric velocity (km s"') 


Fig. 7. Examples of normalized profiles of atomic features as ob¬ 
served at multiple dates in V1309Sco. Most characteristic velocities 
are marked with vertical dashed lines. 


-185 km s'. iMason et al.l (1201 (ill derived the systemic velocity 
of -80+15 km s ' from the position of [Ca ii] lines in their spec¬ 
tra from earlier dates than those considered here. Their value 
corresponds well to the position of emission components in our 
spectra, which - affected by blueshifted absorption - cannot be a 
good measure of the systemic velocity. The corresponding FSR 
velocity of Vsys+SR = -175 km s“' places V1309 Sco at the edge 
of the velocity ra nge seen in the inters tellar emission of H i to¬ 
wards its position (iKalberla et al.l2005h . This would suggest that 
VI309 Sco is a very distant object, perhaps at the outskirts of the 
Galaxy (~20kpc). However, the distance determinations consid¬ 
ered in iTvlenda et al.l (1201 Ih strongly suggest otherwise, giving 
the best estimate of the distance of only 3.5 kpc. From this, we 
conclude that the line profiles in VI30 Sco do not provide the 
true systemic velocity when methods used for classical winds 
are applied. The extreme negative shift may be perhaps related 
to a geometry of the outflow which gives us a good view on the 
approaching parts of the outflow, hiding the redshifted parts. 

The lines do not show any obvious features that would indi¬ 
cate they ar ise in a disk, existenc e of which is evident in other ob¬ 
servations (iTvlenda et al.ll20lTl see also below). For the epochs 
c overed by our data, w e also do not favour the model proposed 
in iMason et al.l (1201 Oh of an expanding shell thicker in the equa¬ 
torial plane and seen from the polar direction. 
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Fig. 5. NIR spectra of VI309 Sco covering the JH bands. Simulated spectra of AlO and VO are shown to indicate the location of the emission 
bands of those species. Location of the water absorption bands is marked with horizontal bars near the edges of the spectral ranges shown. A 
telluric spectrum is also included at the bottom of the figure (with photospheric absorption features of the standard star). Spectra are in arbitrary 

units and are not in scale with respect to each other. _ 
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Fig. 6. The same as Fig.[5]but for spectra covering approximately the HK bands. In addition to AlO, for which a simulation is shown, the first- 
overtone band of CO is clearly present in the 2009 spectra and one rovib line of H 2 is seen in the 2012 spectrum. 
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6. Discussion and conciusions 


6.1. Temporal changes 

The picture of the evolution of VI309 Sco arising from our spec¬ 
tral observations is the following. In 2009, the stellar remnant 
was surrounded by cool molecular gas, but atomic and ionized 
gas, which dominat ed the spectra of V1309 Sco in earlier epochs 
(iMason et al.ll20ldll . was still present. As our simple modelling 
of water bands showed, the molecular gas had a rotational (ki¬ 
netic) temperature lower than about lOOOK, but the gas was 
clearly not in thermal equilibrium as even higher vibrational 
temperatures of up to 3000 K are evident in spectra of other 
species. The considerable strengths of emission features shows 
that the emitti ng area had to be s ignificantly larger than the stel¬ 
lar disk (e.g. ISchoenberall985|) . The molecular material was 
most likely expanding in a form of a shell during whole 2009, 
as can be seen by the absorption features becoming weaker with 
time and by the increasing intensity of emission features (relative 
to the continuum and the corresponding absorption features); a 
wind-like outflow is also evident in the P-Cygni type profiles of 
some atomic lines and molecular features. 

During the whole period covered by our data we clearly ob¬ 
served the B aimer lines and forbidden lines of metals of low 
ionization degree. This is evidence that hot gas of a temperature 
of a few thousand K and low density was present in the system. 
The decreasing flux of these features is due to quick cooling. It 
is most natural to assume that this is the material heated during 
the 2008 eruption and still dissipating its energy. However, the 
observation of the forbidden line of H2 present only in 2012 sug¬ 
gests that the remnant might still have been dynamically active 
as late as 2012 and produced portions of gas of high excitation. 
The emission in [Oi] was optically thin during the whole pe- 
riod covered by our observations, as it was in the earlier epochs 
(iMason et al.ll2010h . 

The source of optical/NIR continuum disappeared sometime 
between the end of 2009 and April 2012 and the remnant of 
VI309 Sco was seen only in discrete spectral emission features 
in our latest observations. The disappearance may be related to 
the decrease in brightness of the stellar remnant owing to its 
cooling but it is likely to be partially related to increased ob¬ 
scuration by dust. 

Dust could already ha ve been present in V1309 Sco in 2008 
just before the eruptio n dTvlenda et al.l 1201 it [hfeCollume^^ 


2014; IZhu et al1l2()13h . but absent before 2008 (iNicholls et ~ 


20131) . As explained in the earlier studies, the dusty material 


formed sometime before the eruption and is likely confined to 
a disk seen almost edge on. The optical emission in continuum 
and lines seen in the 2009 spectra can be interpreted as scattered 
light of a stellar source that had been considerably obscured to 
us by the dusty disk. The scattering would occur on the dust 
(continuum) and associated atomic/mole cular gas (lines), just 
as described in our studies of V4332 Sgr (iKamihski et al.ll2010l: 


iKamiriski & Tvlendall2013t iTvlenda et al.ll2014l) rThis scattering 

effect is even more pronounced in VI309 Sco with a drastically 
low intensity ratio of the Na i to K i doublets. 

The fact that we do not see the continuum source in 2012 
must be related either to the formation of fresh dust that would 
increase the extinction towards the source or to some kind of 
geometrical reconfiguration of the material around the star, for 
instance by developing a thicker flared disk that would block the 
light along our line of sight. Infrared photometry from 2010 sug¬ 
gests the former scenario because it shows a sudden increase of 
redde ning of the source (iMcCollum et al.l 12014t INicholls et al.l 
I2OI3I) . Our spectral observations therefore bracket the (second) 


dust formati on period of VI309 Sco . We note however that, as 
suggested in lMcCollum et al] (l2014l) . the spectral features in the 
NIR spectra are strong compared to continuum and can influ¬ 
ence interpretation of the broadband photometry in the sense that 
the fluxes do not correspond to pure continuum, especially after 
2009. 

6.2. V1309 Sco as a member of the red transient class 

As described in detail in IMason et al.l (l2010l) . the spectral be¬ 
haviour of V1309 Sco before 2009 followed that of another red 
transient, V838Mon, which erupted in 2002. The 2009 opti¬ 
cal spectra presented here still bear strong similarity to those 
of V838 Mon in later stages of its evolution dLvnch et al.ll2004 
HootI), in particular owing to the deep absorption bands of vari¬ 
ous oxides. The infrared bands of AlO , which are rare in astro¬ 
nomical spectra (iBaneriee et al.ll2()T^ but strong in V1309Sco 
in 2009, were also seen in V838Mon, although only in absorp¬ 
tion (Lynch et al.). 

The spectral appearance of VI309 Sco in 2012 is remark¬ 
ably similar to what has been observed in V4332 Sgr years to 
decades after its outburst (1994). In fact, the emission spectrum 
seen at optical wavelengt hs is almost exactly t he same as that 
of V4332 Sgr described in Tylenda et alJ (l2005l) . IKamihski et al.l 
(l2010l) . and iTvlenda et al. (l2014l) . The only difference is the 
presence of the forbidden emission lines of [Oi] and [Feii] 
in VI309 Sco, which, we predict, should have disappeared by 
now (2015). Similarly, we predict that the UVB spectrum of 
VI309 Sco will evolve to what is seen in V4332Sgr, i.e. a spec¬ 
trum dominated by low-excitation lines of neutral atoms and 
molecules. The otherwise very rare infrared features of AlO 
and CrO seen in emission make late VI309 Sco most akin to 
V4332Sgr, as does the silicate dust feature observed near 10/zm 
in both objects (INicholls et al.l201^ and references therein). The 
geometry and excitation conditions in the merger remnants of 
both objects are comparable as well. 

The similarity of late spectra of VI309 Sco to those of the 
other longer-known red transients is striking. The resemblance 
strengthens the mutual link between the older red transients and 
VI309 Sco and suggests similar physical mechanisms driving 
their post-outburst evolution. Because the merger event is best 
seen in VI309 Sco (the Rosetta stone) - owing to a wealth of 
photometric data before and during the spiralling-in process - 
the link provides strong support to the interpretations in which 
all the older Galactic red transients were indeed caused by stellar 
mergers. 


6.3. Chemical composition 

The omnipresence of oxides - in particular the strong water 
bands seen in the 2009 spectra - are evidence that the molec¬ 
ular remnant of VI309 Sco is oxygen rich. This is in agreement 
with the pr esence of inorganic (s ilicate) dust in the remnant as 
observed bv iNicholls et al.l(l2013l) . There are, however, some no¬ 
ticeable chemical/spectroscopic peculiarities in VI309 Sco. One 
is the presence of rare NIR bands of AlO. Their presence in red 
transients a nd a relation to format i on of alumin a dust have been 
discussed in iBaneriee et al.l (l2003ll2004l2()T^ . 

The presence of CrO in VI309 Sco (and V4332 Sgr) is even 
more remarkable. No other circumstellar spectral features of 
CrO have been detected so far, includi ng visual, infrared, an d 
millimetre/sub millimetre wavelengths dSheridan et al.l l2002h . 
The fact that we observe CrO in red transients but not in other 
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known classes of objects implies that either explosive chemistry 
of red transients favours formation of this particular metal ox¬ 
ide or that CrO is exceptionally abundant owing to an increased 
abundance of elemental chromium (or only one of its isotopes, 
see below). Deriving molecular abundances in red transients is 
very challenging and may in fact be impossible until a very good 
knowledge of physical conditions and geometry of their circum- 
stellar surroundings is acquired. Even then, deriving elemental 
abundances from molecular ones would be a very difficult task, 
especially since these eruptive objects are unlikely to be in chem¬ 
ical equilibrium. Although the field of non-equilib rium chem¬ 
istry is starting to achieve its first successes (e.g. ICherchn^ 
l2006l) . much effort is needed before realistic chemical models 
for objects like VI309 Sco will be possible. 

Nevertheless, the unique presence of CrO in red transients 
points to the virtual possibility of abundance anomalies in 
their remnants. An effort should be taken to identify the re¬ 
maining unassigned molecular features in the NIR spectrum 
of VI309 Sco to investigate whether the potential enhancement 
in chromium is associated with overabundance of other rare 
species. The iron-group elements are of particular interest here 
since oxides of other transition metals (TiO, VO, ScO) are 
clearly present in VI309 Sco. Interestingly, it was recently re¬ 
alized that anomalous abundances of ^"^Cr, ^^Fe, and ®^Ni in so¬ 
lar system materials may not originate from supernovae, as ear¬ 
lier assu med, but can be a produ ct of i-processing in low-mass 
stars (e.g. lWasserburg et al.ll2015l and references therein). Some 
collisions of white dwarfs are also expected to produce mate¬ 
rial enhanced in iron-group elements, especially in '^"^Ti and "^^Cr 
(e.g. lPapish & Per^l2015L and references therein). Whether the 
strong signature of CrO in spectra of red novae is a manifes¬ 
tation of a related nucleosynthesis phenomenon is an intrigu¬ 
ing question, especially in the context of the recent discovery of 
rare isotopologues in t he oldest recognized red transient, CK Vul 
(iKamihski et al.ll2015h . 


Appendix A: Simulations of the molecular bands 

For identification purposes and first excitation analysis of the 
observed bands, we performed simple simulations of their elec¬ 
tronic or vibrational bands. More detailed excitation analysis for 
chosen species will be presented elsewhere. All mole cules ex¬ 
cept f or water (see below) were performed in pgopher (IWesternI 
l2014l) . The simulation is very basic and does not include full 
radiative-transfer calculations. In particular, only a very limited 
correction for line saturation effect in the high-opacity regime 
was included. All simulations were performed as for a stationary 
isothermal layer (slab) of gas. Each species was simulated in¬ 
dividually. For some electronic systems, the relative intensities 
of the different rovibrational transitions were calculated using 
Franck-Condon factors, if available in the literature. 

AlO A (NIR) : Spec troscopic constants were taken 

from iFaunila & JonssonI (Il994h (their Table 1) and the simula¬ 
tion included vibrational states u' < 4 in the lower electronic 
state and v” < 8 in the upper electronic state. The Frank-Condon 
factors (FCFs) were taken from Table 3 of [Bernard & Gravinal 
(Il984l) to obtain the relative intensity scale of the bands. Only 
the main isotopologue, ^^Al'^O, was included in the simulation; 
no significant contribution from the unstable A1 isotopes was no¬ 
ticed in the observed spectra. 


VO (VIS): The (visual) B - X system of VO was 

sim ulated with spectrosc opic constants from lAdam et al. I (Il995h 
and ICheung et alJ (Il994h . Those sources provided us with data 
for the 11=0 and 1 in both upper and lower electronic states. 
Some constants were extrapolated to the v''=2 state, allowing 
for a rough identification of the (0,2) band position but result¬ 
ing in high uncertainty in the b and structure. FCFs o f the (0,0) 
and (1,0) bands were taken from lCastelaz et akl (l2000h . For (0,1) 
band, the FCF was assumed to be the same as for (1,0), serving 
only as a very rough approximation on the band intensity. 

VO A^n - and - A^n (NIR): In the NIR region, 
the electronic structure of VO is very complex, in particular be- 
cause of overlaps between bands of different systems (see e.g. 
iMerer et'anil987l) . We simulated the NIR spectrum of VO in¬ 
cluding the (0,0), (0,1), (0,2) bands of the A - X system and the 
(0,0) band of the D-A system. In the laboratory spectrum of VO, 
bands with negative values of At; in the B — X system (described 
in the paragraph above) overlap with the near-infrared systems; 
e.g. the B — X (0,3) band blends with A - X (0,0). We did not 
include those B - X bands from higher vibrational levels in our 
simulation, mainly because they are expected to be weak and be¬ 
cause of a lack of spectroscopic constants for the high vibrational 
levels. For t he bands of th e A-X system, the spectroscopic con¬ 
stants from ICheung et al] (Il982h were used. Our pgopher sim¬ 
ulation of this syste m also included the re lative band intensities 
based on FCFs from lCastelaz et aP (120001) . We also included the 
(0,0) band of the D-A system, which would overlap partially 
with the A-X (0,0) band, but found no indication of its presence 
in the observed spectra. The simulation was made only for the 
main isotopologue, ^'V*®0. 

CrO A^E - X^n (NIR): The A - X (0,0) band of CrO was 
simulated in pgopger using an example file provided by the 
program’s a uthor, C. Western, an d based on spectroscopic con- 
stants from Cheung et al.l (Il984l) . Using spectroscopic data of 
iBarnes et al. (11993 ) for the c = 1 level of the ^11 electronic 
state, we extended the simulation to the (0,1) band, which is pre¬ 
dicted to be located near 1365-1430nm. This places the band 
within the saturated band of water and therefore the observed 
spectra could not be conclusively tested for the presence of this 
CrO band. We also added the (1,0) band by calculating spectro- 
scopic data for the v - 1 state at A ^E from the AG1/2 value from 
ICheung et al.1 (Il984l) . This allowed us to predict the origin of this 
band near 1 lOOnm, a region dominated by absorption of water. 
The position of the strongest head (spike) within the band, how¬ 
ever, overlaps with a discrete emission-like feature seen in all 
2009 spectra, but the shape and the central position are not so 
well matched as in the case of the (0,0) band. We therefore could 
not conclusively decide on its presence. 

CrOA'5A-X5n(VIS): The band A' - X ( 0,0) of CrO was sim- 
ulated using spectroscopic constants from iBarnes et alJ d 19931) 
provided for the v - 0 and 1 states within the upper and lower 
electronic levels. For u = 1 at X^A, the band origins are given 
individually for the five different components arising from the 
spin-orbit interaction. We calculated the corresponding vales 
parametrizing the interact ion (i.e. 7’n,A, A, 77, 0, fixing y at 0.00), 
following formula (19) in ICheung et alJ (Il984l) . With these con¬ 
stants, we performed simulations for the (1,0), (0,0), (1,1), and 
(0,1) bands. However, so far the laboratory spectra have only 
been reported for the overlapping (0,0) and (1,1) bands so our 


Article number, page 12 of[l3] 


































































Kaminski et al.; Stellar merger remnant VI309 Sco 


simulation of the (1,0) and (0,1) bands is rather uncertain. The 
observed shape of the (0,0) band suggests a rotational tempera¬ 
ture of lower than ~300 K and may be as low as a few tens of 
K. 

The predicted position of the (0,1) band near 905 nm over¬ 
laps with the strong pcrr band of water and the (0,2) band of 
VO hampering any conclusive identification. If the vibrational 
temperature is close to the rotational temperature we are able 
to roughly constrain, this band should be very weak. The (1,0) 
should be stronger. The predicted location of the (1,0) band over¬ 
laps with the strong 0,0 band of VO at about 790 nm. There are 
several strong features with no identification within the VO band 
whose shape is similar to that of the simulated (1,0) band of CrO. 
However, the position of the simulated band misses the uniden¬ 
tified emission by ~1.6nm. Better spectroscopic data for vibra- 
tionally excited levels are needed to conclusively decide whether 
the observed emission can be identified as CrO bands. 

H2O (VIS/ NIR): Simu l ations of water bands were based on a 
line list of iBarber et al.l (l2006l) and included only the main iso- 
topologue, H^^O. Unlike for all the other molecules, the simu¬ 
lations of H2O were obtained with the codespectra-BT2 . £90 
program, which is distributed with the line list. The line intensi¬ 
ties are calculated under conditions of thermal equilibrium and 
within the limit of optically thin lines. The spectra were simu¬ 
lated as arising from an isothermal and stationary slab of mate¬ 
rial. 
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